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ABSTRACT 


LINE  CRACK  SUBJECT  TO 
ANTIPLANE  SHEAR1 

A.  Cemal  Er ingen 
Princeton  University 


Field  equations  of  nonlocal  elasticity  are  solved  to  determine 
the  state  of  stress  in  a plate  with  a line  crack  subject  to  a 
constant  antiplane  shear.  Contrary  to  the  classical  elasticity 
solution,  it  is  found  that  no  stress  singularity  is  present  at 
the  crack  tip.  By  equating  the  maximum  shear  stress  that  occurs 
at  the  crack  tip  to  the  shear  stress  that  is  necessary  to  break 
the  atomic  bonds, the  critical  value  of  the  applied  shear  is 
obtained  for  the  initiation  of  fracture.  If  the  concept  of 
the  surface  tension  is  used, one  is  able  to  calculate  the 
cohesive  stress  for  brittle  materials. 


1.  INTRODUCTION 


In  several  previous  papers  [1]  - [4]  we  discussed  the  state  of 
stress  near  the  tip  of  a sharp  line  crack  in  an  elastic  plate  subject  to 
uniform  tension  and  in-plane  shear.  The  field  equations  employed  in  the 
solution  of  these  problems  are  those  of  the  theory  of  nonlocal  elasticity. 
The  solutions  obtained  did  not  contain  any  stress  singularity,  thus 
resolving  a fundamental  problem  that  persisted  over  half  a century.  This 
enabled  us  to  employ  the  maximum-stress  hypothesis  to  deal  with  fracture 
problems  in  a natural  way.  Moreover  it  has  been  possible  to  predict 
the  atomic  cohesive  stresses  by  introducing  the  experimental  values 
of  the  surface  energy. 
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The  present  paper  deals  with  the  problem  of  a line  crack  in  an  elastic 
plate  where  the  crack  surface  is  subject  to  a uniform  anti-plane  shear 
load.  This  problem,  classically,  is  known  as  the  Mode  III  displacement. 

We  employ  the  field  equations  of  nonlocal  elasticity  theory  to  formulate 
and  solve  this  problem.  The  solution,  as  expected,  does  not  contain  the 
stress  singularity  at  the  crack  tip  and  therefore  a fracture  criterion 
based  on  the  maximum  shear  stress  hypothesis  can  be  used  to  obtain  the 
critical  value  of  the  applied  shear  for  which  the  line  crack  begins  to 
become  unstable.  If  the  concept  of  the  surface  energy  is  introduced,  it 
is  possible  to  calculate  the  cohesive  stress  holding  the  atomic  bonds 
together.  Estimates  of  cohesive  stress  are  also  given  for  perfect 
crystalline  solids. 

In  section  2 we  present  a resume  of  basic  equations  of  linear  non- 
local elastic  solids.  In  section  3 the  boundary-value  problem  is  form- 
ulated, and  the  general  solution  is  obtained.  In  section  4 we  give  the 
solution  of  the  dual  integral  equations,  completing  the  solution.  Cal- 
culations for  the  shear  stress  are  carried  out  in  a computer,  and  the 
results  are  discussed  in  section  5. 

2.  BASIC  EQUATIONS  P*1  NONLOCAL  ELASTICITY 

Basic  equations  of  linear,  homogeneous,  isotropic,  nonlocal  elastic 
solids,  with  vanishing  body  and  inertia  forces,  are  (cf.  [5]): 


(2.1) 


tk2, , k 


0 


(2-2)  Hr 


[ A * ( | x * —x  I ) e_„  (x* ) 6,,0+2y  * ( | x * — x | )ek£  (x' ) ]dv  (x' ) 


'rr  C ' k£ 


ekt  ■ i<uk,«+ut,k> 
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where  the  only  difference  from  classical  elasticity  is  in  the  stress 
constitutive  equations  (2.2)  in  which  the  stress  t^(x)  at  a point  x 
depends  on  the  strains  e^^(x’)»  at  all  points  of  the  body.  For  homo- 
geneous and  isotropic  solids  there  exist  only  two  material  moduli, 

A * ( | x * —x | ) and  p'(jx'-x|)  which  are  functions  of  the  distance  |x'-x|. 

The  integral  in  (2.2)  is  over  the  volume  1/  of  the  body  enclosed  within 
a surface  3^. 

Throughout  this  paper  we  employ  cartesian  coordinates  x^  with  the 
usual  convention  that  a free  index  takes  the  values  (1,  2,  3),  and  repeated 
indices  are  summed  over  the  range  (1,  2,  3).  Indices  following  a comma 
represent  partial  differentiation,  e.g. 

V*  ■ 8Va\ 

In  our  previous  work  [6,  7]  we  obtained  the  forms  of  A ' ( | x ' — x | ) 
and  u'(|x-x|)  for  which  the  dispersion  curves  of  plane  elastic  waves 
coincide  with  those  known  in  lattice  dynamics.  Among  several  possible 
curves  the  following  has  been  found  very  useful 

( A ’ , p * ) « (A,y)a(|x-x| ) , 

(2.  A) 

a ( | x * —x | ) = aQexpt-  (3/a)2(x'-x) • (x'-x) ] , 


where  3 is  a constant,  a is  the  lattice  parameter,  and  aQ  is  determined 
by  the  normalization 


a( | x'-x| )dv(x' ) * 1 


(2.5) 


V 


is  the  classical  Hooke’s  law.  Substituting  (2.8)  into  (2.1)  and  using 
Green-Gauss  theorem  we  obtain: 


(2.9) 


ct(|x?-x|)  okJ^k  (x*)  dv  (x')  -1  a ( jx’-x  j)  ok£  (x'Jda^x*)  = 0- 

l'  * 


The  contribution  to  the  surface  Integral  from  the  parts  of  the 
surface  at  infinity  would  be  dropped  since  the  displacement  field 
vanishes'  at  infinity. 


3.  CRACK  UNDER  ANTIPLANE  SHEAR 


We  consider  an  elastic  plate  in  the  (x^«x,  x2«y)  - plane  weakened 

by  a line  crack  of  length  21  along  the  x-axis.  The  plate  is  subjected 

to  a constant  anti-plane  shear  stress  t *t  along  the  surfaces  of  the 

yz  o 

crack.  Fig.  1.  For  this  problem  we  have 
(3.1)  u1-u2«0  , u3-w(x,y)  , 

°-2)  <V"|5  - V'ly  • a11  °ther  V0  ' 

so  that  the  only  surviving  member  of  the  field  equations  (2.9)  is 

(3.3)  yj  a(|x'-x| , |y'-y| )V,2w(x' ,y')dx'dy*-  j o( |x'-x| , |y | ) * [e yz(x' ,0) ]dx'«0 

J l 

where  the  integral  with  a slash  is  over  the  two-dimensional  infinite 
space  excluding  the  line  of  the  crack  (|  x|  <H,  y*C) . A boldface  bracket 
indicates  a jump  at  the  crack  line. 

When  an  undeformed  and  unstressed  body  is  sliced  to  create  a free 
surface,  it  will  in  general  be  deformed  and  stressed  on  account  of  the 
long-range  interatomic  forces.  Thus  if  we  are  to  consider  that  the 
Plate  with  a crack  is  undeformed  and  unstressed  in  its  natural  state 
then  we  must  apply  the  boundary  conditions  on  the  unopened  crack  surface. 

Under  the  applied  anti-plane  shear  load  on  the  unopened  surfaces 
of  the  crack, the  displacement  field  possesses  the  following  symmetry 
regulations 

(3.4)  w(x,-y)  - -v(x,y,  . 


Using  this  in  (3.2)  we  find  that 


(3.5)  [o  (x,0)]  ■ 0 
- y* 

Hence  the  line  integral  in  (3.3)  vanishes.  By  taking  the  Fourier 
transform  of  (3.3)  with  respect  to  x* , we  can  show  that  the  general 
solution  of  (3.3)  is  identical  to  that  of 

(3.6)  - C2w(S,y)  ■ 0 , 

dyZ 

almost  everywhere.  Here  a superposed  bar  indicates  the  Fourier 
transform  e.g. 

f(£.y)  - (2ir)“^  J"f(x,y)  exp  (iCx)  dx  . 

The  boundary  conditions  are 

v(x,0)  « 0 for  j x | >1  , 

0*7)  ^z**'0)  " T0  for  |x|  <£. 

w(x,y)  - 0 as  (x2+y2)1/2^»  . 

The  general  solution  of  (3.6)  (for  y>0)  satisfying  (3,7)3  is 

(3.8)  v(x,y)  - (2/ir)^r  A(0  e”£y  cos  «x)  d£  , 

'o 

where  A(£)  is  to  be  determined  from  the  remaining  two  boundary  conditions. 
For  the  'non-zero  components  of  the  stress  tensor  we  have 
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(3.9) 


Cxz  “ -(2/w>i  V £AU)«e  j"dy»  Pi  o ( | x ' -x j > I y * -y | ) 
-o( jx'-x| , |y'ty| )]  e”^y  Sin  (£x')  dx’ 5 
tyz  - -<2/*)4  y | A(0*de  pdy*  pi  a(|x’-x|,iy’-y|) 
+a(|x'-x| , (y'+y|)]  e“^y  cos  (£x')  dx'. 


Using  (2.4>2for  a(jx'-x|),  we  carry  out  integrations  on  x'  and  y'.  To 
this  end  we  note  the  following  integrals,  18]: 

i, . [■ «p(-p«'2>  j | *•  - <’/P>*  «xp(-62/aP)  (sgj  \ , 

(3.10)  I2  - pexp(-py'2-Yy')dy'  • ±(7r/p)*  exp(Y2/4p) [l-« (Y /2^p)  , 

'o 

i: 


#(z)  = 2ir  3 1 exp(-t2)  dt  . 
' o 


Hence 


(3.11)  t 


y* 


- -(270"*  y 

r 5A«)  [e"5y  erf cl 

^~2py\  _ e^y  erfc 

* 

0 

V 2/?  / 

V 2/p  / 

- -(2ir)"*  y j 

P £A(£)  I«”5y  erfc 

(&n\  + ,«:■  ,rfc 

/5i2£zV 

1 0 

\2/p) 

\ 2/?  / 

(8/a)2 

erfc(t)  - 1- 

♦ (*)  . 

The  boundary  conditions  (3.7)^  and  (3.7)2  now  read 


z-x/fc,  C ■ Ci  , e ■ a/231  . 
(3-13)  K(eC)  ■ erfc  (eC)  » 

A(0  - C_*  C(C)  , Tc  - Xq12/v  . 


I 

i 

I 

\ 

\ 


\ 

i: 


To  determine  the  unknown  function  A(0 » we  must  solve  the  dual  Integral 
equations  (3,12). 


THE  SOLUTION  OF  THE  DUAL  INTEGRAL  EQUATIONS 


Recalling  the  e*preszion 


r I 

I cos  (zc)  » (itz;/2 )*  (z;)  , 

P 

i 

[ where  J^z)  is  the  Bessel  function  of  order  v,  we  write  the  system 

| (3.12)  in  the  form 


The  kernel  function  k(c()  Is  given  by 


(4.2)  k(e;)  « K(e;)  -1  - -*(c;)  • 

The  solution  of  the  dual  Integral  equations  (4.1)  is  not  known.  However, 
it  is  possible  to  reduce  the  problem  to  the  solution  of  a Fredholm 
equation  (cf.  [9]) 

1 

(4.3)  h(x)  + f h(u)L(x,u)du  - -i(itx)^TQ  * 

o 

for  the  function  h(x),  where 

(4.4)  L(x,u)  - (xu)*|  tk(ct)Jo(xt)J  (ut)dt  • 

'o 

Vhen  (4.3)  is  solved,  then  C(;)  is  calculated  by 

(4.5)  0(0  - (20*  [ x*JQ(Cx)  h(x)  dx  • 

'o 

As  discussed  in  a previous  work  [4],  if  we  note  that  e is  extremely  small, 

k(eO  ®*y  be  neglected  as  compared  to  unity  In  (4.1),  (see  Fig.  2),  In 

this  case  the  zeroth  order  solution  of  (4.3)  namely  h (x)  • -T  (nx)*/2 

o o 

suffices  for  the  calculations  vhen  the  crsck  size  is  larger  than  100  atomic 
distances.  In  such  a case  we  have 

CQ(0  - -(ir/2)^To;“*J1(;) 


(4.6) 
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and  therefore 

(4.7)  Ao(0  - -(ir/2)iToJ1(5£)/a 

The  shear  stresses  are  then  calculated  by  (3.11).  Interesting  among 

these  is  the  shear  stress  t along  the  crack  line  y » 0.  For  this 

yz  J 

we  obtain 

(4.8)  t (z,0)/to  - j K(eO  J1(;)  cos(zO  dc 

y J 0 

As  observed  before,  this  integral  converges  for  all  z provided  K(e£) 

is  not  approximated  by  unity  for  e small.  For  e * 0 at  z ■ 1 we  have 

the  classical  stress  singularity.  However,  so  long  as  e 4 0,  (4.8) 

gives  a finite  stress  all  along  y * 0.  At  0<z<l,  tyz/T0  is  very  close 

to  unity,  and  for  z>l,  t /tq  possesses  finite  values  diminishing  from  a 

maximum  value  at  z=l  to  zero  at  z=°° . 

For  e>>!/100  the  approximate  solution  given  by  (4.7)  is  not  very 

good.  However,  further  improvements  can  be  achieved  by  the  iterative 

solution  of  (4.3)  with  the  use  of  Cq(C).  Since  e> 1/100  represents  a 

crack  length  of  less  than  10  ^cm,  and  at  such  submicroscop ic  sizes  other 

serious  questions  arise  regarding  the  Interatomic  arrangements  and  force 

laws,  we  do  not  pursue  solutions  valid  at  such  small  crack  sizes. 

5.  NUMERICAL  CALCULATIONS  AND  DISCUSSION. 

Calculations  of  the  shear  stress  t , given  by  (4.8)  along  the  crack 

yz 

line,  were  carried  out  on  a computer.  The  results  are  plotted  for  e **  1/20, 
V 50,  1/100,  1/200,  in  Figures  3 to  6.  For  a crack  length  of  20  atomic 


distances  (e  - 1/20)  the  result  is  not  very  good  in  that  the  boundary 


-11- 


condition  at  jxjcl,  y ■ 0 is  satisfied  only  very  roughly.  However,  for  a 
crack  size  of  100  atonic  distances  (Fig.  5)  the  shear  stress  boundary 
condition  io  fulfilled  in  a strong  approximate  sense.  The  relative  error 
in  this  case  is  less  than  li°/o.  Hence  we  conclude  that  the  classical 
Ao(£)  given  by  (4.7)  gives  satisfactory  results  for  crack  lengths 
greater  than  100  atomic  distances. 

The  stress  concentration  occurs  at  the  crack  tip,  and  this  is  given  by 

(5.1)  t U,0)/t  = c /Je  , e - a/201 

y*  o 3 

where  c^  converges  to  about 

(5.2)  c3-  0.40  . 

The  following  observations  are  very  significant; 

(i)  The  maximum  shear  stress  occurs  at  the  crack  tip,  and  it  is 
finite  (eq.  5.1) 

(il)  The  shear  stress  at  the  crack  tip  becomes  infinite  as  the  atomic 
distance  a + 0.  This  is  the  classical  continuum  limit  of  square 
root  singularity. 

(iii)  When  t ^(1,0)  ■ t£  ( ■ cohesive  shear  stress) , the  plate  will  fail 


In  this  case 
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Equation  (5.3)  is  the  expression  of  the  Griffith  fracture  criterion 
for  brittle  fracture.  We  have  arrived  at  this  result  via  the  maximum 
shear-stress  hypothesis,  rather  than  the  surface  energy  concept  used 
by  Griffith  and  his  followers.  The  significance  of  this  result  is  that 
the  fracture  criteria  are  unified  at  both  the  macroscopic  and  the 
microscopic  scales  and  that  the  natural  concept  of  bond  failure  is 
employed . 

(iv)  The  cohesive  shear  stress  t^  may  be  estimated  if  one  employs 
the  Griffith's  definition  of  the  surface  energy  y and  writes 

(5.5)  t 2a  - K y > 
c c 

where 


r P 

j;  (5.6)  Kc  «=  8yc3e/v(l-v) 

i- 

r 

; Since  some  measurements  exist  on  y , by  employing  these  values  we 

I 

jf  can  calculate  the  cohesive  shear  sttess.  For  steel  we  have 


I 

ft.' 


r 

r 


r 

E 


Y = 1975  CGS  , w = 6.92  x 1011  CGS 

v = 0.291  , a = 2.1+8  A0 

(5.7)  tc/p  = 0.2568  B1/2 
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Let  at  a n atomic  distance  the  nonlocal  effects  attenuate  to  1$  of  its 
value  at  x=0.  Using  (2.1+)  we  find  that 

(5.8)  3 = 2.lH6/n  , 

(5.9)  t /y  = 0.3761+  n_1/2  . 

c 

For  n = 6 this  gives 

(5.10)  t /y  = 0.14 

c 

which  is  in  the  right  range  and  well  accepted  by  metallurgists.  For 

example  Kelly  [10]  gives  t /y  = 0.11.  There  is  however  a question  of 

c 

in-plane  versus  antiplane  shear  failure  which  need  special  examination. 
Compared  to  the  results  obtained  in  our  previous  work  on  the  in-plane 
shear  [1+],  the  cohesive  stress  seems  30%  higher.  However,  this  is  somewhat 
artificial  since  it  is  necessary  to  know  the  value  of  3 or  n in  either 
case  more  precisely.  This  of  course  requires  at  least  one  experiment. 

Acknowledgement : The  author  is  indebted  to  Dr.  Balta  for  carrying  out 
computations. 
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